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THE ZERO -LIFT WAVE DRAG OF A FAMILY OF OPEN-NOSED 
BODIES OF REVOLUTION IN THE MACH NUMBER RANGE 
OF 2.0 TO 4.0 

By Leroy L. Presley and Emmet A. Mossman 
SUMMARY 


The wave drag of a family of open-nosed bodies of revolution was 
computed by six approximate theories, first-order perturbation theory, 
second-order perturbation theory, generalized shock-expansion theory, 
second-order shock-expansion theory, tangent-wedge theory, and impact 
theory, and by the method of characteristics for Mach numbers between 
2.0 and 4.0. The best agreement with the method of characteristics was 
provided by second-order perturbation theory and second-order shock- 
expansion theory with the latter being the most attractive from the 
standpoint of computing time required versus accuracy obtained. 

The wave drag, for bodies of the type investigated in this study 
with the same initial lip angle and the same diameter ratio (initial to 
maximum diameter) , was found to increase sharply for fineness ratios less 
than 3 but to be nearly constant for fineness ratios above 3. The wave 
drag was found to decrease nearly linearly as the diameter ratio increases 
for bodies having the same initial lip angle and fineness ratio. 

An approach is given for adapting the method of characteristics to 
automatic computing machine procedure. 


INTRODUCTION 


For aircraft configurations employing a pod-nacelle arrangement to 
house an air breathing propulsion system, the external wave drag of the 
engine housing can be a significant portion of the total drag of the 
aircraft. In the absence of systematic experimental data to aid in 
estimating this wave drag, the usual approach is to calculate it theoret- 
ically. At present there are a number of different theoretical methods 
which can be used for such calculations. However, the results obtained 
from the various methods are not consistent over a wide range of Mach 
numbers and body shapes and the computation time differs widely among 
the methods. The investigation described herein was made, therefore, 
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to compare the theoretical wave drag as computed using six approximate 
methods with that from the method of characteristics in order that a 
choice from the viewpoint of average computing time required and 
accuracy obtained could be made. The theories which were compared to 
the method of characteristics were: first-order perturbation theory 

(refs. 1 and 2), second-order perturbation theory (ref. 2), generalized 
shock -expansion theory (refs. 3 and 4), second-order shock-expansion 
theory (ref. 5), tangent-wedge theory (ref. 6 ), and impact theory (ref. 7) • 

To provide a reference for the comparison, the study was conducted 
for a family of open-nosed bodies of revolution having a fineness ratio 
of 5, a ratio of initial to maximum diameter of 0.742, and lip angles 
between 1.478° and 35.844°. This family was selected since its contours 
are representative of the nacelles of present day pod-nacelle arrange- 
ments for jet engines. Thus the wave-drag characteristics of such a 
family are of interest in themselves. The effects of varying fineness 
ratio and diameter ratio were also determined using the method of charac- 
teristics and generalized shock-expansion theory. 


THEORETICAL METHODS 


Although the various theoretical methods used in the present study 
have been discussed in detail in their respective references, it is con- 
sidered appropriate to this report to discuss briefly each of the methods, 
pointing out in particular their approach to the solution of the super- 
sonic flow field about an open-nosed body of revolution and their expected 
range of applicability. All of the theories used in this analysis have 
as their basis a solution of the gasdynamics equation shown below: 



bv _ uv /chi 
Sr a 2 \dr 




( 1 ) 


where the symbols are defined in appendix A. Since this equation, which 
is applicable to any steady inviscid flow of a perfect gas, is nonlinear, 
simplified methods of solution must be used for most problems. The 
simplified methods of solution which have resulted in the theories used 
in this study are: numerical solution of equation (l), linearization 

of equation (l) with subsequent analytical solution of the linearized 
equation, and approximate solutions of equation (l) which are applicable 
to certain flow regimes. In the application of the theories, the follow- 
ing conditions were imposed: 

(1) The flow entered the nose of the body at supersonic speed (i.e., 
mass-flow ratio of unity). 

(2) The bodies were immersed in an ideal gas. 

(3) The bodies were at zero angle of attack. 
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Method, of Characteristics 


One method of solving the gasdynamics equations is hy numerical 
means. Such a method is greatly facilitated in the present case of 
supersonic flow since the gasdynamics equation is a hyperbolic differ- 
ential equation and of a type that is integrable on characteristic sur- 
faces which correspond to Mach lines. It is then necessary to rewrite 
equation (l) in a form suitable for numerical integration along the Mach 
lines, and several different forms have been derived (refs. 8, 9 , and 10). 
One form of the compatability equation for the method of characteristics 
which was considered to be the most satisfactory for the present study 
since it contains only two flow- field variables is as follows: 1 


with the upper sign referring to the first family Mach line and the lower 
to the second family Mach line as shown in figure 1. Equation (2) 
expresses the relationship between static pressure, stream angle, and 
Mach angle along Mach lines in the flow field and is applicable to either 
rotational or irrotational steady flow of an isoenergetic gas. When 
applied to rotational flow, as in this study, the change in entropy 
normal to the streamlines is taken into account by considering the change 
in total pressure along the Mach lines. This method is limited to bodies 
with supersonic flow behind the nose shock. 

In the actual solution of the flow by the method of characteristics, 
the compatability equation is put into finite difference form and solved 
point for point in the flow field. If solved by purely numerical means, 
the computations are very long and laborious for any practical problems 
and become feasible only with the use of automatic computing machines. 

In the present study, the computations were made using an automatic 
computing machine. The equations used as well as a discussion of some 
considerations in the adaption of the method of characteristics to 
automatic computing machine procedure are given in appendix B. 


In contrast to the numerical integration of the nonlinear gasdynamics 
equation by the method of characteristics, first-order theory introduces 
iThis form is not given directly in either references 8, % or 10 
but can be obtained by a suitable combination of equations 3*21 to 3*23 
in Chapter I of reference 10. 



( 2 ) 


where 


a = -hE 


sin 2 [i 


(3) 


First-Order Theory 


k 
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a perturbation function into the equation and linearizes the resulting 
expression so that it may be solved analytically. This is done by con- 
sidering the magnitude of the velocities in the following manner 


u = U + Up 
Up and Vp « U 


v = v T 


u-n and v-n « a 




and defining a perturbation potential such that 


ox 


Vp dr 


Equation (l) can be put into the following form 


1 - U + 


a 2 cp 

dx 2 


1 - 

\dr/ 


d 2 cp + 1 d£ _ 2 d^ dep d 2 cp 
dr 2 r dr dx dr dx dr 


( 5 ) 


If the squares and products of derivatives of the perturbation potential 
are neglected, the following equation results 


dr 2 dx 2 r dr 


( 6 ) 


where (3 2 = M 2 -l. Equation (6) is the linearized form of the potential 
flow equation. Details of the integration of this equation are given in 
references 1 and 2. Such a solution is applicable to the calculation of 
flow fields at values of the hypersonic similarity parameter, Moa& 0 , less 
than 1.0. It in effect neglects entropy losses due to the presence of 
the body since entropy losses are of third order in the perturbation 
potential. 


Second-Order Theory 


The third approach to obtaining a solution of equation (l) is to 
use an iteration procedure, an approach first considered by Busemann 
and later extended by Van Dyke (ref. 2). For this method, a first-order 
solution is obtained as described previously. This solution is substi- 
tuted into the right-hand side of equation (5) and a second-order 
perturbation potential is found. The method of solution is given in 
reference 11. Therein it was stated that the initial angularity of the 
body must be less than the free-stream Mach angle. 

In the present study the computations were done using an automatic 
computing machine since they are lengthy when done by hand. The initial 
angularity of the bodies was restricted to I3 0 or less because of the 
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limitations in the automatic computing machine program, an angle con- 
siderably less than that permitted by the theory. This theory is also 
applied only in the range of Mx>So < 1 . 0 . 


Generalized Shock -Expansion Method 


A method which has as its basis an approximate solution of the exact 
equation of motion is the well-known generalized shock-expansion theory. 

The approximation made for this theory consists of neglecting the term 
(ds/r)sin p sin S of equation (2) for cases in which the rate of change 
of surface angle of the body is large compared to its divergence. The 
resulting expression is seen to be a differential form of the Prandtl -Meyer 
equation. The flow field downstream of the nose shock wave can then be 
considered of the Prandtl-Meyer type. The aforementioned assumption 
regarding the body shape is made in this theory so it could be expected 
to become accurate when MoJS>o » 1.0 and when the area ratio of the body 
(maximum to initial area) is near 1 . 

The application of this method begins by first approximating the 
body by a series of straight-line elements, tangent to the original body. 
The flow at the nose is defined by means of the Rankine-Hugoniot relations 
for an oblique shock wave. Since the flow downstream of the nose is of 
the Prandtl-Meyer type, the pressure can be found on any tangent line if 
the pressure and Mach number on the preceding tangent-line element are 
known. The pressure is taken to be constant on any given tangent-line 
element. The method is limited as is the method of characteristics to 
bodies that have supersonic flow behind the nose shock wave. The calcu- 
lations for this method can be made efficiently using a desk calculator 
and the tables and charts of reference 12 . 

In the present investigation, the body contour was approximated by 
13 straight-line elements. 


Second-Order Shock -Expansion Theory 


Second-order shock-expansion theory was developed to provide a 
closer approximation to the flow field for the cases when M^Sq is near 
1.0. This was done by developing an approximation to the part of equa- 
tion (2) that was neglected for the generalized method. Two significant 
differences from the generalized method arise out of this closer approxi- 
mation. First, the exact pressure gradient at the nose of the body is 
found from the method of characteristics. Secondly, along each tangent 
line used to approximate the body contour, the pressure is found to vary 
exponentially. The asymptote of the exponential pressure variation is 
assumed to be equal to the pressure on a cone having the same slope as 
the tangent- line element and at the same free-stream Mach number. 
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The application of this method is similar to that of the generalized 
method except the pressure gradient on each tangent line is found hy the % 

relations given in reference 5. The results can he obtained by means of 
a desk calculator but the process is somewhat tedious. 

In the present investigation, the body contour was approximated by 
13 straight-line elements. 

Tangent -Wedge Method 

An approximation to the generalized shock-expansion method is the 
tangent-wedge method. This approximation is based on the assumption that 
the pressure on the body at any point is a function of the net flow deflec- 
tion angle from the free stream. Thus in this method, the pressure at any 
point on the body can be found from the Rankine-Hugoniot relations for 
oblique shock waves (for bodies inclined into the free stream) or the 
Prandtl-Meyer relations (for body surfaces inclined away from the free 
stream) in which the deflection angle and the upstream Mach number used 
in these relations are the local angle of the body and the free-stream 
Mach number, respectively. The method is therefore limited to bodies with 
angularity less than the detachment angle of a two-dimensional shock wave 
and could be expected to become accurate in the same regions as the gener- 
alized shock-expansion method. The computations for this method are very * 

simple and can be made efficiently by means of a desk calculator and the 
charts of reference 12. 

In the present investigation, the pressure was computed at Ik points 
along the body. 


Impact Theory 

A direct approximation to the tangent-wedge method for certain flow 
fields has resulted in defining an area where Newtonian impact theory 
becomes applicable. For flows at infinite Mach number, wherein 7 = 1.0 
and the shock wave is coincident with the body contour, the expression 
for the pressure coefficient at any point on the body as given by the 
tangent-wedge method can be shown to reduce to 

Cp = 2 sin 2 5 

This expression was obtained by Newton by neglecting centrifugal forces and 
assuming that the component of the momentum of the free-stream air that is 
normal to the body surface is absorbed, thereby creating a force on the 
surface. This theory which is applicable to any body contour would be 
expected to become applicable for M 00 ^ ) » 1.0. Computations for this theory 
are very simple and can be made efficiently using a desk calculator. 

In the present investigation, the pressure was computed at Ik points 
along the body. 
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COMPUTATIONS 
Body Contours 


Througtiout this investigation, a family of profiles whose contours 
are a function of initial lip angle, S 0 , the diameter ratio d^/dm and 
the fineness ratio l/dm was used. The equation relating the bodies is 
given in reference 13 and has the following form: 

r^ = r m - (rm-roHl-x) 1 " 1 (8) 

where r\ is the relating parameter for the family of curves and is 
defined hy 


T = 


tan 5 0 



( 9 ) 


The diameter ratio, c^/c^, was varied from 0.707 to O. 898 . However, 
for the majority of the investigation, a value of 0.7^2 was used. This 
diameter ratio is in the range of values considered in reference It and 
of those necessary to envelop existing turbojet and ramjet engines at 
about Moo = 3*°* Th e diameter ratio was varied while the initial lip 
angle and fineness ratio were held constant. 


The majority of the theoretical investigation was conducted with 
profiles having a fineness ratio of 5* These profiles are shown in fig- 
ure 2. A small part of the theoretical investigation was concerned with 
varying the fineness ratio of the bodies from 0 .625 to 10 while the 
initial lip angle and diameter ratio were held constant. 


Wave -Drag Computation 


The wave drag was obtained from the following relationship 

C D =/0pd (10) 

The methods discussed previously were used to calculate Cp as a function 
of A^/Aq. The above integration was carried out graphically. 

RESULTS AND DISCUSSION 
Pressure Distributions 

The distribution of the external pressure coefficients as computed 
by the seven theories described earlier is presented in figures 3 through 
9, and are in the form of static-pressure coefficient, Cp, as a function 
of local area ratio, A^/Aq. These curves were subsequently integrated 
to obtain the wave-drag coefficients. 
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Comparison of typical pressure distributions is shown in figure 10. 

From these curves three general comparisons can he made as to the agree- 
ment between the method of characteristics and the various approximate 
theories. The first point of comparison is the magnitude of the initial 
pressure rise at the nose of the body. Generalized shock-expansion, 
second-order shock-expansion, and tangent-wedge theory (omitted from the 
figures because of its closeness to generalized shock-expansion theory) 
have the same initial static -pres sure coefficient as the method of 
characteristics since all of these theories use the exact pressure at 
the nose. First-order perturbation theory and impact theory give lower 
values of initial static -pres sure coefficient than the method of charac- 
teristics whereas second-order perturbation theory gives a higher value. 

The second point of comparison is that only second-order perturbation 
theory and second-order shock-expansion theory predict an initial pres- 
sure gradient which is similar to that predicted by the method of 
characteristics. It should be mentioned that the initial pressure 
gradient of the second-order shock-expansion method is inherently iden- 
tical to that of the method of characteristics. Finally, the predicted 
variation of static -pres sure coefficient downstream of the nose differs 
for the various theories. For a curved body such as investigated herein, 
the method of characteristics predicts positive pressure coefficients 
over most of the body and an overexpansion with resulting negative static - 
pressure coefficients near the base of the body. Three of the theories, > 

first- and second-order perturbation theories, and second-order shock- 
expansion theory, compare favorably with the method of characteristics 
for predicting negative static -pres sure coefficients for curved bodies. ▼ 

However, positive pressure coefficients were predicted over the entire 
body by generalized shock-expansion theory, tangent-wedge theory, and 
impact theory. To summarize , the results of figure 10 indicate that 
second-order perturbation theory, and second-order shock-expansion theory 
provide the best agreement with the method of characteristics for these 
body shapes and Mach numbers. 

The variation of the static -pres sure coefficient for bodies of 
varying fineness ratio (i/dm of O.625 to 10), but with constant initial 
lip angle and diameter ratio, is shown in figure 11 for = 2.5. These 
distributions were computed by the generalized shock-expansion theory 
(fig. 11(a)) and the method of characteristics (fig. 11(b)). Both 
theories indicate that at a given area ratio (A^/Aq), a lower static 
pressure can be obtained with a body of higher fineness ratio. The 
significance of this will be discussed later. 

The effect upon the static-pressure coefficient distributions of 
varying the diameter ratio, while the initial lip angle, fineness ratio, 
and M^, are held constant, is shown in figure 12. It can be seen that 
the pressure distributions for the bodies are similar in that at the same 
proportionate area, the pressure is approximately the same. 


NACA TN 4368 


9 


WAVE DRAG 


The theoretical wave-drag coefficients for the bodies having a 
fineness ratio of 5 and a diameter ratio of 0.742 at Mach numbers of 
2.0, 2.5, 3 . 0 , and 4.0 are presented in figure 13 . As indicated, there 
is a large difference in the results as obtained from the various methods 
throughout the Mach number range of the investigation. These differences 
are seen better in figure 14 in which the wave drag computed by the method 
of characteristics is used as a reference and the error in the wave drag 
as computed by the six approximate methods is shown. The error is given 
as a function of the two-dimensional hypersonic similarity parameter, 

M» 6 o. It can be seen that both generalized shock-expansion theory and 
tangent-wedge theory overestimate the wave drag while impact theory 
underestimates the wave drag throughout the entire range of Moo 5 0 
investigated. The error for each of these methods is sizable for values 
of <1.0 but decreases as M^&q increases. The two perturbation 

theories have good accuracy at low values of M^d 0 (l^ 8 0 = 0.1 to 0 . 2 ) 
but start becoming inaccurate as M^o increases. The iteration in 
second-order perturbation theory is seen to be effective in increasing 
the range of accuracy over that of the first-order solution. Second- 
order shock-expansion theory has good accuracy in the range near 

= 1 but becomes inaccurate at extremely low values of M^b. The 
apparent regions of applicability of the various theories are as would be 
expected from their assumptions discussed previously. A summary of the 
wave drags obtained by the method of characteristics is presented in 
figure 15 in the form of a wave-drag parameter, C D (O.7M00 2 ), as a function 
of Moo & o (ref. 15) • It is seen that the wave drags correlate well with 
a curve representing the mean of the values at any given M^Sq. 

The effect of varying the fineness ratio, l/d TO , on the wave drag as 
computed by the method of characteristics and generalized shock-expansion 
theory^is shown in figure 16 for bodies with an initial lip angle of 
21.156 and a diameter ratio of 0.742, and at a Mach number of 2 . 5 . 

As was mentioned earlier, the bodies with higher fineness ratios had a 
lower static pressure existing at the same area. Thus these bodies would 
have a lower wave drag. The results as shown indicate that the wave drag 
of open-nosed bodies of revolution (mass flows of unity) is a function of 
the fineness ratio and that a fineness ratio of at least 3 is necessary 
to obtain near minimum wave drag. 

The variation in wave drag with changes in diameter ratio is shown 
in figure 17 for generalized shock-expansion theory and the method of 
characteristics. The data show that the difference in wave drag as 
computed by the two theories decreases as the diameter ratio approaches 
1 , a result that could be expected from the assumptions of generalized 
shock-expansion theory. The wave drag is also seen to decrease almost 
linearly as the diameter ratio increases. 
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Computing Time 


The final choice of which theory to use in any given case must he 
governed by the accuracy desired and the computing time available. The 
computing time and accuracy of the various theories are shown in the 


following table. 

Theory 

Method of characteristics 
First-order perturbation 
Second-order perturbation 
Generalized shock expansion 
Second-order shock expansion 
Tangent wedge 
Impact 


Computing time. 


hr 

Accuracy 

40 

Best 

16 

Fair 

24 

Good 

2 

Poor 

8 

Good 

1 

Poor 

1 

Poor 


The computing times given are for one person using a desk calculator. 

All solutions except the method of characteristics are for 14 points on 
the body. The computing time given for the method of characteristics is 
for a solution in which the mesh size was allowed to become large 
(c/r = 1.0) . Such a solution will give wave-drag coefficients slightly 
less than those presented in this study. If greater accuracy is desired 
(finer mesh size) the computing time for hand computation can be considered 
to increase in proportion to the square of the ratio of the mesh size for 
c/r = 1.0 divided by the mesh size vised (c/r < 1.0). 


CONCLUDING REMARKS 


► 


> 


The wave drag of a family of open-nosed bodies of revolution was 
computed by six approximate theories and the method of characteristics 
for a Mach number range of 2.0 to 4.0. Using the results from the method 
of characteristics as a reference, the investigation showed that the 
three theories which required the shortest computing time, generalized 
shock-expansion, tangent -we dge , and impact theories were in general the 
least accurate. First- and second-order perturbation and second-order 
shock-expansion theories gave more accurate results, predicting the wave 
drag to within 10 percent of that from the method of characteristics over 
much of the ranges of Mach number and body shapes investigated. The range 
of applicability of the theories was found to vary with the perturbation 
theories giving good results at Moo5 0 < 1.0, second-order shock expansion 
giving good results near Moo5 0 1.0 with the generalized shock-expansion, 

tangent -we dge, and impact theories starting to give good results at 
Moo&o » 1.0. In general, considering computing time as well as accuracy, 
second-order shock-expansion theory was the most attractive within the 
range of variables investigated herein. 
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The effect of varying the fineness ratio , 2 /c^, on the wave drag was 
investigated "by the use of two theories, the method of characteristics 
and generalized shock-expansion theory. The results of doth theories 
indicate that for bodies of the type investigated in this study and having 
the same initial lip angle and diameter ratio, near minimum wave drag can 
be obtained with a fineness ratio above 3 . 0 . 

The wave drag was also computed for bodies having various diameter 
ratios with the same initial lip angle and fineness ratio. For these 
bodies, both the method of characteristics and generalized shock-expansion 
theory predict a near linear decrease in wave drag as the diameter ratio 
increases. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Aug. 21, 1958 
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APPENDIX A 


NOTATION 


For convenience the list of symbols has been divided into two groups, 
those used throughout the paper and those used exclusively in appendix B. 
Following are those symbols in the first group. 


A 

a 

C D 


D 

d 

l 

M 

P 

Pt 

q 

r 

s 

U 

u 

v 

X 


area, sq ft 

speed of sound, ft/sec 

wave- drag coefficient, — i- 

OdoA. 

static-pressure coefficient. 


Pb-P 0 

Qoo 


distance from point to point along Mach lines in method of 
characteristics solution divided by l 

drag, lb 

body diameter, ft 
total length of body, ft 
Mach number 

static pressure, lb/sq ft 
stagnation pressure, lb/sq ft 
dynamic pressure, lb/sq ft 
radial distance divided by l 

distance along Mach lines in method of characteristics 
free- stream velocity parallel to x axis 
axial component of velocity 
radial component of velocity 

distance along axis of body from origin divided by l 
n/m 2 - 1 

ratio of specific heat at constant pressure to specific heat 
at constant volume 
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5 stream angle, radians 

T) parameter relating body shapes (see eq. (9)) 

6 shock- wave angle, radians 

X -&K- 

sin 2|i 

\i Mach angle, sin 1 i, radians 

M 

p mass density, slugs/cu ft 

cp perturbation potential 

Subscripts 

b body 

m maximum geometric characteristic of body 

o characteristics of body at origin 

p perturbation quantities 

00 free-stream conditions 


The following symbols are used exclusively in appendix B which lists 
the equations used in calculating the wave drag by the method of character- 
istics. 


Cp distance along first family Mach line from computed point to 

immediately preceding upstream point divided by l 

c s distance along second family Mach line from computed point to 

immediately preceding upstream point divided by l 

n number of the computed point 

(The number sequence proceeds from point nearest body to shock 
wave and from ray to ray in downstream direction (see fig. 1).) 

N number of points in input ray (N = 5 in fig. 1) 

e a small arbitrary number indicating closeness of approximation 

in iterative methods 
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f 


n 

s 


w 


Subscripts 

characteristics at point whose number is n-1 and 'which 
lies on first family Mach line immediately preceding 
computed point 

characteristics at point of number n 

characteristics at point whose number is n+l-N and which 
lies on second family Mach line immediately preceding 
computed point 

characteristics at point whose number is n-N and which 
lies on shock wave immediately preceding computed point 
on shock wave 


i, 2 ,...i number of iterations 


Superscripts 


iterations 
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APPENDIX B 

CALCULATIONS BY METHOD OF CHARACTERISTICS 


Solutions by the method of characteristics presented in this paper 
were computed by an automatic digital computing machine. The equations 
used in the solution are also applicable to manual computation. They 
are therefore listed and discussed herein for ready reference of the 
interested reader. 

In the application of the method of characteristics to the calcu- 
lation of the static pressure on a body of revolution, the axially 
symmetric flow field bounded by the body surface and the nose shock wave 
is subdivided by a network of lines inclined to the local streamline 
at the local Mach angle (see fig. 1). Points are defined as the inter- 
section of a pair of such lines (field points), the intersection of a 
line with the body contour (body point), or the intersection of a line 
with the shock wave (shock-wave point). A ray is defined as the line 
containing a series of points connected by first family Mach lines. 

Ifru-S in figure 1, points 1 to 'J and 8 to 14 lie on two separate rays. 

Two types of equations are required in the solution, those which 
define the coordinates of the point and those which define the aero- 
dynamic properties at the point. The calculation of the coordinates 
and aerodynamic properties proceeds from point to point along a ray, 
commencing at the body and ending at the shock wave, and then from 
ray to ray in the downstream direction until the end of the body is 
reached. In the following discussion, the equations will be grouped 
according to whether the point is on the input ray, or a field, body, 
or shock— wave point downstream of the input ray. The equations are 
given in terms of an arbitrary point so they can be used in the 
repetitive type of calculation required by the method of characteristics. 


Characteristics of Points on Input Ray 


input ray is located on a first family Mach line sufficiently 
close to the nose of the body that the flow properties at the points can 
be obtained with the assumption that the shock wave is two-dimensional. 
In the present case the input ray intersected the body at x b < 0.001. 
The coordinates of points on the input ray are - 


x, 


•n 


N-n 

N-l 


x b 


n-l r^-*btan(p + 5)- r Q- 
N-l _ tan 0-tan(n+5) ^ 


(Bl) 
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and 

_ N-n n-1 f[r) 3 -Xbtaxi(p.+&) ]tan 0-r o tan(|i+&) 

n N-1 ^ N-1 \ tan 0-tan(|i+&) 

■where 9 is the two-dimensional shock-wave angle corresponding to a 
deflection equal to S Q and p and 6 are the flow characteristics 
behind the wave. 

It should be mentioned that the equations discussed hereinafter are 
adaptable for determining the entire flow field about a closed-nosed 
body of revolution as well as the open -nosed bodies studied in this paper. 
In the former case, the flow behind the nose shock must be supersonic. 



Characteristics of Field Points 


The equations for the coordinates of the field points are: 

Xo-Xf +r s cot ( p-6 ) s+rf cot ( p+S ) f 

r n = — — 

co t ( |i- & ) s +co t ( p+5 ) f 


(B3) 


and 


*n = (r s -r n )cot(p-5) s +x s 


(B4) 


The distances along the Mach lines from the calculation point to the 
preceding known points are given by: 


c f = 


c s - 


rg-rf 

sin(p+5)f 

r s -rn 


(B5) 


sin(p-5) s 

The first and second family compatibility equations, equation (2) 
in the body of the paper, are put in finite difference form expressing 
the difference between the aerodynamic properties at point n and the 
preceding points along the first and second family Mach lines, and are 
solved simultaneously to give for the stream angle at point n 


(B6) 


&n = 







- tf- + A5 



AS - 
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and the static -pres sure ratio 



c s (sin p sin 6) 


0 


(b8) 


Assuming the total pressure to vary linearly with distance normal to 
streamline and the variation of stream angle to he small in the vicinity 
of point n gives 



n 


opsin \if 


CfSin [i f +c s sin p. s 




Jr (B9) 


The ratio of local static pressure to local total pressure is 


from which 


and hence 





. -1 

sin 


JL_ 

M n 


Characteristics of Body Points 


(BIO) 


(Bll) 


(B12) 


The body point n is at the intersection of the second family Mach 
line from point n + 1 - N and the body contour. Hence the coordinate 
equations are dependent on the analytic expression for the body contour. 
If the egression is not simple, an explicit solution of the intersection 
point may be impossible, as in the present case. In such cases Newton’s 
approximation can be vised to solve for the axial location of the point 
as follows. Let the difference between the body radius and the radial 
distance of the second family Mach line from point n + 1 - N at the 
same value of x be expressed as 

f(x) = r b (x)-r s +(x-x s )tan(p-6) s (B13) 

where rj-,(x) is the analytic expression for the body contour, so that 
at point of intersection 


(B14) 


V 


f(xn) = 0 
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To solve for let x = X! = x s and solve 


x 2 = x - 


f(x) 


= Xl - 


f( x i) 


df (x)/dx * dfCx-LJ/dx 

Repeat this operation letting x 2 x 2 and continue until 

xi-xi_i < e 


(B15) 


(B16) 


The value of x n is taken to be xp. Then the radius at the inter- 
section point is 

rn = r s -(xn-x s )tan(p-S)s (B17) 

The distance along the second family Mach line from point n to 
the preceding point, c s , is given by equation (b 6) and the stream angle, 
5 n , is given by the slope of the body contour. These two values are used 
in equation (B8) to determine the static pressure on the body at point n. 
The static -pres sure coefficient is then given by 


, = (P/Pjn-B 

p (r/2)Moo 2 


(B18) 


Since the body contour is a streamline, the total pressure remains 
constant along the body. The remaining flow characteristics are 
determined from equations (BIO) to (B12). 


> 




Characteristics of Shock-Wave Point 


The shock-wave point, n, is at the intersection of the first family 
Mach line from point n-1 and the shock wave from point n-N. The 
coordinates are given by 

x w tan 0 w -Xftan(|i+5)f+rf-r w 

x n ; 

tan 0 w -tan(p+&)f 

and 

r n = ( x n" x f ) tan ( p+5 ) f +rp 

The distance along the first family Mach line, Cf, is given by equa- 
tion (B5) . 

A simple explicit relation involving the static pressure and 
deflection angle immediately behind a shock wave is not available. It 
is therefore impractical, and probably impossible, to determine the static 


(B19) 


(B20) 
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pressure and deflection angle at the shock-wave point by solving simul- 
taneously the equation relating pressure and deflection angle across the 
shock and the one compatibility equation available. Iterative methods 
are therefore used. The compatibility equation available is 


5 n 


V) JjL) 

Pojf Wn «, cf (sin p sin S)f 

— + of - • ■ 

Af r f 


To commence the iterative procedure, let 



(B21) 


(B22) 


and find 6 n ' from equation (B21) . The stream angle 5 n ' is taken 
as the deflection angle for a two-dimensional shock wave at point n 
and the shock-wave angle computed from 


sin 6 0 n '+g sin 4 0 n '+h sin 2 0 n '+j = 0 


(B23) 


where 


. 7 sinV 

Moo 


h = 


2MJ+L 


(7+ 1 ) 2 . (7-1) 


Mx^ 


sin 2 ^ 1 




J = - 


cos 2 ^* 

M» 4 


(B24) 


Equation (B23) is solved most readily by Newton's approximation. 
Briefly let 

f(d') = sin 6 0’+g sin 4 0'+h sin 2 0'+j (B25) 

so that at point n 

f(e') = f(0 n ') = 0 (B26) 

To solve for 0 n ' , let 0’ =03/ = 0 W and solve 

sin20 2 * = sin 2 0 * sin 2 0 1 ' f (V ) 

df (0 ' )/d sin 2 0 1 dfCei’J/d sin 0 ! ’ 

(B27) 
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Repeat the operation letting 0 ' = d 2 ' and continue until 

Qi’-^i-l' < e (B28) 

The value of 0 n ' is taken to he 9±' . 

The static-pressure ratio across a shock wave corresponding to a 
given shock- wave angle is 


_p 27Hx> 2 sin 2 0-(7~l) 

Poo 7+1 

The static pressure at point n is then taken as 


(B29) 


Let 



1 /jp \ ' + 27H3o 2 sin 2 0-(7-l) 

2 \Poo/n + 7+1 


(B30) 



n 


P_ 

Poo/n 


11 


(B31) 


in equation (B21) and solve for 6 n " . The process involving equations 
(B23) to (B31) is repeated until 

Sn^-^n^ -1 < e (B32) 


Then 



(B33) 


The total pressure loss through the shock wave is given by 

7 

7-1 




(7+1)Mqo sin 0 ^ 
(7-l)H» 2 sin 2 0n+2 


P 

Pocr 


1 

1-7 


(B34) 


The remainder of the aerodynamic characteristics are found by means of 
equations (BIO) to (B12) . 
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Mesh Size 


Accurate solutions of the method of characteristics require that 
the distance between points, mesh size, be small. However, when the 
mesh size, c/r, becomes large, it is necessary to iterate the points 
affected. This iteration consists of averaging the aerodynamic properties 
of the known and computed points, assigning them to the coordinates of the 
known point and recomputing the computed point. In the actual machine 
computation, a c/r of 0.25 was used as the criterion for iteration when 
N = 10 and 0.10 when N = 20, and the iteration was performed only once. 

The computing time (on an IBM 653) for 10 and 20 point solutions 
was 3/4 and 3 hours, respectively. For curved bodies, the difference in 
static -pressure coefficient distribution from 10 and 20 point solutions 
was insignificant. For straight bodies, however, it was found necessary 
to start at x < 0.001 in order to obtain a pressure distribution with 
no discontinuities. 
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Figure 1.- Schematic diagram of method of characteristics solution. 
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Figure 2.- External contours of fineness -ratio -5 bodies. 
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(a) Hx, = 2.0 

Figure 3*~ The variation of static-pressure coefficient with local-area 
ratio; method of characteristics; fineness ratio = 5.0 ; diameter 
ratio = 0.7^2. 
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(b) Moo = 2.5 

Figure 3- - Continued. 
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(c) Mpg = 3-0 
Figure 3 • - 


Continued 
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(d) = 4.0 

Figure 3 .- Concluded. 
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Figure 4.- The variation of static -pres sure coefficient with local-area ratio; first-order 
perturbation theory; fineness ratio = 5*0) diameter ratio = 0.742. 
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Static-pressure coefficient, 
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Figure 5.- The variation of static -pres sure coefficient with local-area ratio; second-order 
perturbation theory; fineness ratio = 5.0; diameter ratio = 0.742. 
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(a) M^, = 2.0 

Figure 6.- The variation of static -pres sure coefficient with local-area 
ratio; generalized shock-expansion theory; fineness ratio = 5.0; 
diameter ratio = 0.7^2. 


coefficient , C 


34 


NACA TN 4368 



Loco I- area ratio , /A 0 


(Id) = 2.5 

Figure 6.- Continued 
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(c) = 3.0 

Figure 6.- Continued 
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Local-area ratio , A b /A 0 

(d) Mqq = ^.0 

Figure 6.- Concluded. 
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(a) M = 2.0 

Figure 7 .- The variation of static-pressure coefficient with local-area 
ratio; second-order shock-expansion theory; fineness ratio = 5 .O; 
diameter ratio = 0.742. 
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(b) M = 2.5 
Figure 7 .- Continued 
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(c) M = 3.0 


Figure 7 .- Continued 



Static pressure coefficient, Cp 
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(d) M = 4.0 

Figure 7* - Concluded. 
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(a) = 2.0 

Figure 8.- The variation of static -pressure coefficient with local-area 
ratio; tangent-wedge theory; fineness ratio = 5.0; diameter 
ratio = 0.742. 
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(b) Moo = 2.5 
Figure 8.- Continued 
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(c) Moo = 3.0 
Figure 8.- Continued 
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(d) Moo = 4.0 

Figure 8.- Concluded 
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Figure 9 .- The variation of static-pressure coefficient with local-area 
ratio; impact theory; fineness ratio = 5 * 0 ; diameter ratio = 0 . 7 ^ 2 . 




46 


NACA TN 4368 



1.0 1.2 1.4 1.6 1.8 2.0 

Local-area ratio , Ajj/Aq 

(a) = 2.0, n = 9 

Figure 10.- Comparisons of typical static-pressure-coefficient distribu- 
tions for bodies having a fineness ratio of 5-0 and a diameter ratio 
of 0.742. 
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(8) H* = 2.0, T) = 15 
Figure 10.- Continued. 
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1.0 1.2 1.4 1.6 1.8 2.0 

Local-area ratio , A b /A 0 

(c) Moo = 2.5, T) = 9 

Figure 10.- Continued. 
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Figure 10.- Continued. 
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Local-area ratio, At/A 0 

(e) ^ = 3-0, t) = 18 
Figure 10.- Continued. 


Static-pressure coefficients 


NACA TN 4368 


51 



(f) ^ = 4.0, ti = 18 


Figure 10.- Concluded. 


Static- pressure coefficient , C 



(a) Generalized shock-expansion theory. 


(b) Method of characteristics. 


Figure 11.- The variation of static -pres sure coefficient with local-area ratio for bodies of 
varying fineness ratio; 5 0 = 21.156°; Moo = 2.5; diameter ratio = 0.7^2. 
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Wave-drag coefficient, C 
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(a) Hx, = 2.0 

Figure 13.- The variation of the external wave-drag coefficient as a function of the lip-angle 

parameters; fineness ratio = 5-0* 
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Figure 13 .- Continued, 
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Wave-drag coefficient, 





(c) Hx, = 3.0 


Figure 13*- Continued 
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Figure 13 .- Concluded. 
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Error in wave drag 
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Figure 14. - Error of the various approximate theories with the method of characteristics as a 
reference for bodies having a fineness ratio of 5-0 and a diameter ratio of 0.742. 
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Figure 15 .- Correlation of wave-drag coefficients obtained by the method of characteristics 

for dQ/dm = 0.742 and l / dm = 5 . 0 . 
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Figure 16.- The variation of external wave-drag coefficient as a function of 
fineness ratio for 5 0 = 21.156°; ^ = 2.5; diameter ratio = 0.742. 
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Figure 17.- The variation of external wave-drag coefficient as a function of 
diameter ratio for 5 0 = 21. 156° and M = 2.5i fineness ratio = 5.0. 
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